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Abstract 

Bose-Einstein correlations (BEC) between final state particles in the reaction 
e"'"e~ — >■ W"'"W~ — > qiq2q3q4 have been studied. Data corresponding to a total 
integrated luminosity of 550 pb~^, recorded by the DELPHI detector at centre- 
of-mass energies ranging from 189 to 209 GeV, were analysed. An indication 
for inter-W BEC between like-sign particles has been found at the level of 2.4 
standard deviations of the combined statistical and systematic uncertainties. 
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1 Introduction 

Correlations between final-state particles in high energy collisions have been exten- 
sively studied during the last decades. They can be due to phase space, energy-momentum 
conservation, resonance production, hadronisation mechanisms, or be dynamical in na- 
ture. 

In the particular case of identical bosons the correlations are enhanced by the Bose- 
Einstein effect [1,2]. These Bose-Einstein correlations (BEC) are a consequence of quan- 
tum statistics. The net result is that multiplets of identical bosons are produced with 
smaller energy-momentum differences than non-identical ones. 

Several aspects of BEC have been measured in hadronic Z decays and are well under- 
stood [3] . It is natural to expect the same behaviour in the hadronic decay of a single W. 
It is, however, not clear how BEC manifest themselves in a system of two hadronically 
decaying W's, in particular between bosons coming from different W's (inter- W BEC). 

The separation between two W's before their decay is of the order of 0.1 fm, compared 
to a typical hadronisation scale of several fm. Therefore, due to the large overlap between 
the two hadron sources, inter-W BEC cannot be a priori excluded. However, it is unclear 
whether these are of the same type as BEC measured inside a single decaying W, where 
they are, in contrast to the traditional Hanbury-Brown and Twiss [4] picture, not related 
to the total hadronisation volume. 

Together with colour reconnection [5,6], the poor understanding of the inter-W BEC 
effect introduces a large systematic uncertainty in the measurement of the W mass in 
the channel e+e^ W"'"W~ qiqiqaOl [7,8]. The current statistical uncertainty of the 
combined LEP measurement in this channel amounts to 35 MeV [9] , to be compared with 
the total systematic uncertainty in this channel of 107 MeV, which is, however, expected 
to decrease with improved measurements of colour reconnection. The effect of possible 
inter-W BEC amounts to 35 MeV [9]. It is thus clear that a better understanding of the 
phenomenon would help in reducing this uncertainty. 

Measuring inter-W BEC is challenging in practice because of a low sensitivity to 
the effect. This is mainly due to the small fraction of relevant particle pairs coming 
from different W's. Moreover, its isolation from BEC inside a single W requires careful 
attention and needs to be as model-independent as possible. 

The scope of this paper is the model-independent analysis of the correlations of like- 
sign hadron pairs in e+e~ W+W~, where both W's decay into hadrons, with the aim 
of determining the presence and size of inter-W BEC. 

The outline of the paper is as follows: In section 2 the mathematical formalism ap- 
plied throughout this analysis is specified and a brief overview of the analysis is given. 
In section 3 experimental details, such as the detector setup and WW event selection, 
are presented. Section 4 focuses on the mixing procedure employed in order to construct 
an inter-W BEC-free reference sample from events where one W decays leptonically. 
Section 5 clarifies the details of the Monte Carlo models for comparison to the data. 
In section 6 a detailed overview of the numerical analysis of the measured correlation 
functions is given including a construction of weights applied in order to increase the 
sensitivity of the analysis, the subtraction of background and the determination of statis- 
tical errors and correlations of the bins. Moreover the parametrisation of the correlation 
function is discussed. In section 7 results are presented and in section 8 the systematic 
uncertainties are discussed. Finally, sections 9 and 10 discuss the results and conclusions 
are given. 
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2 Analysis method 

The mathematical method used to extract a possible inter-W BEC signal is largely 
based on [10] and [11]. In the case of two stochastically independent hadronically decaying 
W's, the single and two-particle inclusive densities obey the following relations: 

p^^il)=p^-il)+p^-il), (1) 

pWW(l, 2) = (1, 2) + (1,2)+ (l)pW" (2) + pW^ (2)pW" (1), (2) 

where p^(l) denotes the inclusive single particle density of one W andpW(l,2) the inclu- 
sive two-particle density of one W. The densities p^^(l) and p^^(l, 2) then correspond 
to the single and two-particle inclusive densities of a fuUy-hadronic WW event. Assuming 
that the densities for the W"*" and the W~ are the same, which is correct if one does not 
look at the absolute sign of the particles' charges, equation (2) can be re-written as 

p^'W(i^ 2) = 2pW(l, 2) + 2pW(l)pW(2). (3) 

The terms p^^(l, 2) and p^(l, 2) can be measured in fuUy-hadronic and semi-leptonic 
WW decays respectively. 

A pair or two-particle density p^^(l,2) is trivial to construct. The correlation mea- 
surement is made difficult by the fact that only the correlations between particles coming 
from different W's are of interest and there is no way of determining where the particles 
originated from. Finally, in order to obtain a correlation function, it is necessary to con- 
struct a reference sample of events without BEC between particles coming from different 
W bosons. This sample corresponds, in our case, to the product of the single particle 
densities pW(l)pW(2). 

Events where only one of the W's decays hadronically can be used to address these 
challenges. Taking two of these independent hadronically decaying W's and mixing them 
to form one event allows an emulation of a fuUy-hadronic WW event, having BEC inside 
each of the W's. By construction these events will have no correlations in pairs from 
different W's and the measurement becomes a direct comparison between two event 
samples, without any model dependence. The event mixing should follow closely the 
electroweak production of WW events. Possible biases of the mixing procedure can be 
estimated by applying the same procedure to large samples of simulated events. Hence, 
the term p^(l)p^(2) in equation (3) replaced by a two-particle density pj^^, obtained 
by combining particles from two hadronic W decays taken from different semi-leptonic 
events. The details of this "mix ing" procedure are explained in section 4. Expressed in 
the variable Q = \/—{pi — ^2)^, where pi,2 stands for the four-momentum of particles 1 
and 2, equation (3) can be re- written as 

p^(g) = 2pW(g) + 2pWW(g). (4) 

Keeping in mind that equation (2) was formulated for independent W decays, test 
observables can be constructed to search for deviations from this assumption. Such devi- 
ations will indicate that particles from different W decays do correlate. The observables 
considered are: 

Ap(g) = pww(g) - 2p^^'(g) - 2p^^{Q), (5) 

^^^^ = 2pW(g)+2pWW(g)- (6) 
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Given the definition of tlie genuine inter- W correlation function Si{Q) [11], it can be 
shown that 



Ap(g) 



If no inter-W correlations exist, the variables Ap{Q) and Sj{Q) will be zero for all 
values of Q, while D{Q) will be equal to one. Inter-W BEG will lead to an excess at 
small values of Q. 

The selection of particles and pairs is straightforward, with the strongest requirement 

that they should originate from the primary interaction. Moreover, the selected WW can- 
didates have a significant background which must be subtracted using a model dependent 
procedure. 

BEG in Z decays have been extensively measured and constitute a natural basis to 
compare with inter-W BEG. The correlation functions measured in Z events use simulated 
events without BEG as reference samples. They are therefore close to being genuine 
correlation functions but with large model-dependent systematic errors and some dilution 
due to particles which are either not pions or which do no originate from the primary 
interaction. When the inter-W correlations are measured it is natural to compare to the 
Z and single-W data using the same fitting functions. Since the inter-W measurement 
uses data as reference, the model dependence is no longer present. 

The mixing procedure, which allows events to be mixed more than once, leads to 
a rather involved description of the statistical properties of the correlation function. 
However, the same mixing can be used to investigate the sensitivity to the inter-W BEG 
effect. The applied mixing reuses semi-leptonic events up to 20 times, which affects the 
precision depending on whether pairs are constructed by mixing or come from inside single 
W's. Finally, a pair- weighting technique was devised which improved the sensitivity and 
is described in section 6.1. For this purpose the mixed reference sample was used to 
determine statistically whether particles come from the same or different W's. 



3 Experimental details 

3.1 The DELPHI detector 

The DELPHI detector configuration for the LEP2 running evolved compared to the 
one at LEPl [12,13]. The main changes relevant to the analysis described in this paper 
were the extension of both the vertex and the inner detectors. This ensured a very good 
track quality also in the forward region down to small polar angles. During the operation 
of the detector in the latter part of the year 2000 one sector of the TPG malfunctioned 
and the data from this period are excluded from the results. 

In order to verify that a track originates from the primary interaction it was required 
that the TPG participated in the measurement of the track. This effectively required the 
track to be within the polar angle region 20° < 9 < 160°. The reconstructed charged 
particles were required to fulfill the following criteria on the momentum, p, the momentum 
error, Ap/p, and the impact parameters with respect to the event vertex in the plane 
transverse to the beam, e±, or parallel to the beam, ey: 

• 0.2 GeV <p< Pbeam; 

• Ap/p < 1; 

• e_L < 0.4 cm; 
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Year 


1998 


1999 


2000 


(GeV) 
C (pb-i) 


189 
158.0 


192 
25.9 


196 
76.9 


200 
84.3 


202 
41.1 


204-209 
163.4 



Table 1: The integrated luminosities, £, for the various years of LEP2 data-taking, 
expressed in units of pb^^. The corresponding centre-of-mass energies are also given. 

• e\\ < 1.0 cm/sin^. 

The two track reconstruction efficiency in DELPHI drops for opening angles below 
2.5°. Since the mixing procedure does not necessarily reproduce this drop in efficiency 
all particle pairs having an opening angle below 2.5° were omitted in all two-particle 
density distributions. These requirements lead to a typical efficiency of about 85% and 
reduce the total fraction of secondary tracks to about 5%. Secondary tracks are typically 
tracks from secondary decays {K^, A°, etc.) or from secondary interactions in the beam 
pipe and with detector material. Particles not coming from the primary interaction or 
not being pions will dilute the observed correlation. The combined effect was estimated 
to reduce the measured BEG to about 70% of the nominal one. This dilution was not 
corrected for due to model dependence and affects all pair densities in nearly the same 
way. When results from different experiments are combined it will be necessary to apply 
such corrections in order to get comparable results. 

3.2 Selection of WW events 

The total analysed dataset amounts to an integrated luminosity of 549.6 pb~^, collected 
with the DELPHI detector during the years 1998-2000. A summary of the integrated 
luminosity per energy point is given in table 1. 

The samples of fuUy-hadronic and semi-leptonic events required for the WW BEG 
analysis were selected using neural networks, developed in [14] and [15]. 

For the fully-hadronic event selection, it was demanded that the events fulfill the 
following requirements: a large enough charge multiplicity, a large effective centre-of-mass 
energy, large visible energy and four or more jets. 

The final selection was performed using a neural network trained on thirteen event 
variables. The dominant background contribution comes from the qq(7) events. All 
other backgrounds are negligible. Hadronically decaying ZZ events, which constitute 5% 
of the selected sample, were treated as signal as they, except for events where at least one 
Z decays into b-quarks, will have similar space-time kinematics. A comparison between 
data and simulated events of the neural network output for the fully-hadronic selection 
is shown in figure 1. 

By requiring a neural network output larger than a given value, a desired purity or 
efficiency can be reached. The whole analysis was repeated for several cuts on the neural 
network output, selecting samples with an increasing purity, ranging from 83% to 97%. 
This allowed the choice of an optimal working point, minimising the sum of the statistical 
and background uncertainty, corresponding to a selection efficiency and purity of 63% and 
92% respectively, with 3252 events selected in total. 

The semi-leptonic events were selected by requiring two hadronic jets, a well-isolated 
identified muon or electron or (for tau candidates) a well-isolated particle associated with 
missing momentum possibly from the neutrino. The missing momentum direction was 
required to point away from the beam pipe. Dedicated neural network trainings were 
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Figure 1: The neural network output variable for the fully- hadronic event selection. The 
light shaded histogram are the signal, while the dark shaded histograms correspond to 
the background processes. The optimised selection cut is indicated by the arrow. 



used for all lepton flavours. Combining all three lepton flavours, an overall efficiency and 
purity of respectively 58% and 96% was reached, corresponding to 2567 selected events. 
The three neural network outputs, corresponding to the three lepton flavours are shown 
in figure 2 for data and simulated events. 

The WPHACT [16] generator with the JETSET [17] hadronisation model was used for 
the simulation of all signal and four-fermion background events. The qq(7) background 
was simulated using the KK2F [18] generator and also hadronized with JETSET. 



4 Mixing procedure 

The mixed two-particle density, pj^^, was constructed by combining the hadronically- 
decaying W's from pairs of different semi-leptonic WW events, qqlrY, from which the 
lepton was removed and irrespective of the charge of the W's. 

The momentum of each hadronic W can be constructed as the visible W momentum or 
the W momentum after a constrained fit imposing energy- and momentum-conservation 
and constraining the two W masses in the event to be equal to 80.35 GeV. After mixing 
the W's, the first method gives mixed events which have smaller missing momenta than 
the fully-hadronic events while the second gives larger missing momenta. It was therefore 
decided to use the average of the visible W momentum with a weight of 0.4 and the fitted 
W momentum with a weight of 0.6 to obtain the final W momentum. This procedure 
gave the best agreement with respect to the missing momenta, and it was cross checked 
that the mixing quality does not depend significantly on the used weights. 

In WW events the W's are nearly back-to-back due to momentum conservation. In 
mixed events this was accomplished by requiring that one W polar angle lay within 10° 
opposite to the polar angle of the other W. Pairings where the two W polar angles were 
within 10° after inverting the z-component of one W were also accepted. The momenta of 
the W's were then approximatively balanced by rotating one W around the beam axis so 
that the W's became back-to-back in the plane transverse to the beam axis. The above 
transformations refiect the azimuthal and forward-backward symmetry of the DELPHI 
tracking detectors. 
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Figure 2: The neural network outputs for the three semi-leptonic event selections. The 
muon channel (a) and electron channel (b) have small background contaminations, due to 
the clear identification of the isolated lepton. (c) The taus are more difficult to identify, 
resulting in a higher background rate. All signal events are shown by the light shaded his- 
tograms, the background events correspond to the dark shaded histograms. The selection 
cuts are indicated by the arrows. 

All mixed events were subjected to the same event selection as the fuUy-hadronic 
events. The agreement between fully-hadronic events and mixed events was verified for 
several event variables and single-particle distributions. Small differences in the distri- 
butions are taken into account in the estimation of the systematic errors. Examples 
are shown in figure 3, comparing the particle momenta, charge multiplicity, total miss- 
ing momentum and reconstructed W mass between simulated fully-hadronic events and 
simulated mixed events. These events were generated with the BEI model described in 
section 5. 



5 Monte Carlo models 

All Monte Carlo generated events were hadronized using the JETSET algorithm 
unless stated otherwise. BEC were included using the local reweighting algorithm 
LUBOEI [19,7]. It takes as starting point the hadrons produced by the string frag- 
mentation in JETSET, where no Bose- Einstein effects are present. Then the momenta 
of identical bosons are shifted such that the inclusive distribution of the relative sepa- 
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Figure 3: Comparison between simulated fully-hadronic events and simulated mixed 
events for (a) the charged particle momenta, (b) the charge multiplicity, (c) the missing 
momentum and (d) the reconstructed W mass. In all plots the WPHACT generator using 
the BEI model was used. 



ration Q of identical pairs is enhanced by a factor f2{Q) > 1, parametrised with the 
phenomenological form 

UQ) = 1 + Xexp{-Q'R''), (8) 

where Q is the difference in four-momentum of the pair, A and R are free parameters 
related to the correlation strength and the spatial scale of the source of the correlations. 

The corresponding change in the momentum of the particles is not unique. In addition, 
energy and momentum cannot be simultaneously conserved. In the model, the momentum 
is always conserved and afterwards all three-momenta are rescaled by a constant factor, 
close to unity, in order to restore energy conservation. Even when BEC is only allowed 
for pairs coming from the same W (BEI), this global rescaling introduces unreasonable 
negative shifts in the reconstructed W mass. The BE32 variant of LUBOEI overcomes 
this problem by including extra momentum shifts to restore total energy conservation, 
instead of a rescaling of the momenta. 

The BE32 model was tuned to hadronic Z decays, keeping all fragmentation parameters 
fixed, giving a satisfactory result for all hadronic Z events and hadronic Z events with 
reduced b-content. The resulting LUBOEI parameters for the correlation strength. A, 
and the correlation length scale, R, are PARJ(92)=1.35 and PARJ(93)=0.34 GeV"^ (= 
0.6 fm), respectively. Monte Carlo sets exceeding ten times the size of the WW data set 
were simulated at each centre-of-mass energy. 
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Figure 4: The two-particle correlation function for a single W decay, obtained from 
semi-leptonic WW events in data and the BE32 Monte Carlo model. The model was 
tuned using Z data. Fits with equation (11) are superimposed. 



The tuned model gives a good description of DELPHI'S Z data (see figure 10(a)) 
and the hadronic decay of single W's. The latter is illustrated in figure 4, where the 
two-particle correlation function defined as 

P iVjMCnoBE 

is shown for selected scmi-lcptonic W decays in data and MC simulation. The dip around 
Q = 0.5 — 1.0 GcV in the BE32 curve in figure 4 is understood to come from the conserva- 
tion of the total multiplicity in the model and is taken into account by the fit in equation 
(11). The signal at low Q values is naturally compensated by a depletion at higher Q 
values. 

In this paper the label BEl is used for the LUBOEl model including only BEC between 
particles from the same W, BEA is used when all particles are subjected to BEC and 
BEG when all BEC are switched off. 



6 Numerical analysis 

The numerical analysis of the BEC measurement is complicated by statistical correla- 
tions introduced by the methodology. Each charged particle occurs in several pairs and 
each semi-leptonic event is used to produce several different mixed events. All these statis- 
tical correlations were included in the numerical analysis and the performance evaluated 
using resampling techniques. 

6.1 Pair weights 

The sensitivity of any inter-W BEC measurement is hmited by the small fraction of 
particle pairs coming from different W's, resulting in a small Q value. This is illustrated 
in figure 5(a), where the fraction of pairs from different W's, denoted as F[Q), is shown. 
It drops to around 15% at very low Q values. 
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Q [GeV] p(Q,y,e ) 

Figure 5: a):The fraction of pairs coming from different W^s, F{Q), obtained for a BEI 
MC sample witliout pair weiglits (dashed line) and using pair weights (full line). b):The 
estimated pair purity, p{Q, 7, 9*), obtained from mixed semi-leptonic events. 

In addition to Q, the Lorentz factor, 7, and the decay angle of a particle pair, 9*, are 
sensitive to whether two particles come from the same or from different W's. The decay 
angle, 9*, is defined as the angle between the momentum vector of one of the particles 
in the two-particle rest frame and the vector sum of the two particle momenta in the lab 
frame. As such, each individual pair of tracks can be estimated to have a probability 
p{Q,'y,9*) to come from different W's. p{Q,'J,9*) was parametrised for this analysis 
using large samples of simulated mixed semi-leptonic W decays. Figure 5(b) illustrates 
the distribution of p{Q,^,9*) for pairs with Q < 0.5 GcV. The BEI model is shown for 
the mixed and same W's and compared to the data results. 

A particle is combined with the other particles in one event when constructing p^{Q). 
It is combined with the other particles in many other events when constructing p^^{Q). 
Therefore, p^^{Q) has smaller local fluctuations than p^{Q) even though they are con- 
structed using the same particles. Using a detailed error analysis it was determined that 
pairs from the same W contribute a factor 1 + c more to the final variance of the BEG 
measurements at low Q-values than pairs from different W's. For this analysis a value of 
c = 1.9 was determined and used in the following. 

The contribution to the statistical variance was estimated for samples of pairs with a 
given purity and was found to be proportional to 1 + c ■ (1 — p(Q, 7, 9*)). Finally, all pairs 
were weighted by p{Q,j,9*) divided by the above variance factor. Using these weights, 
the improvement in the statistical error on the final measurement was 9%. This is the 
reason for choosing the analysis which weights all particle pairs with their information 
content for the final result. The above procedure not only reduces the statistical error 
but makes the analysis more sensitive to pairs which originate from different W's. 

The two-particle densities in Q for the combined data set are shown in figure 6 for 
both like-sign particle pairs and unlike-sign particle pairs, with and without pair weights. 
In both the fuUy-hadronic and semi-leptonic samples, the number of unlike-sign pairs is 
higher than the number of like-sign pairs at Q values below 2 GeV. This is due to the large 
number of resonance decays with masses in this range. The region around Q = 0.7 GeV 
is dominated by tt'^tt' pairs coming from the p resonance which is abundantly present 
in hadronic decays of the W. Reflections of three-body decays are also present in the 
like-sign distributions. The two-particle densities of like-sign and unlike-sign pairs for 
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Figure 6: The two-particle densities, p^^{Q), P^{Q) and Pmi^iQ)^ ^oi like-sign and 
unlike-sign particle pairs, with no pair weights applied (a-c), and including pair weights 
(d-f). The background contribution from Z(7) — ^ qq(7) events in the p^^{Q) distribu- 
tion is also shown. 



the mixed events coincide, the reason for this being that all pairs in this distribution 
contain particles from different events. When pair weights are applied, contributions 
from resonance decays are down-weighted, therefore the like-sign and unlike-sign spectra 
become more similar. 



6.2 Background subtraction 

The histograms in figure 6 show the contribution from qq(7) background events as they 
are simulated with the BE32 model. The density p^^{Q) is, consequently, corrected for 
this background using the expression 

where iVtot and A^^qq are the total number of selected events and the number of selected 
background events, respectively, and ntot and nqq the respective number of particle pairs 
from these events. 

The correlation strength parameter, PARJ(92), was re-tuned to a value of 0.9 in order 
to get a better description of four-jet Z events, which are more similar to the selected 
background than inclusive Z decays. This re-tuning was only used for the background 
events. Details on the background tuning are discussed in section 8. 



6.3 Fit parametrisation 
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sample/parameter 




R{fm) 




<5jv 


xVndf 


BE32 


0.77 ± 0.02 


0.59 ± 0.01 


-0.78 ± 0.02 


0.013 ± 0.003 


141.8/96 


data 


0.64 ± 0.07 


0.59 


-0.78 


0.018 ± 0.017 


114.7/98 



Table 2: Results of the fit, using equation (11), to R2{Q) — 1, obtained from semi-leptonic 
WW decays, both for data and the BE32 model. 



In order to quantify an excess at small Q values, fits were performed to the inter- 
W correlation function 6i{Q). The choice of fitting function is motivated by the shape 
of the i?2(Q) measurements. Si{Q) and i?2(Q) have nearly the same physical meaning 
but different systematics since i?2(Q) is dependent on the fragmentation model used. 
This means that the optimal fitting functions are not necessarily the same and that 
comparisons between R2{Q) and Si{Q) results must be done carefully. 

It is known from BEC measurements of the hadronic final state of a Z or a single W that 
two particle correlation functions are reasonably well described by either a Gaussian or an 
exponential parametrisation [3,20]. However, the BE32 Monte Carlo samples show a more 
detailed structure in the Q range between 0.5 and 1.5 GeV, as can be seen from figure 4. 
The slope observed around 1 GeV flattens out above Q = 2 GeV. Therefore in the plots 
of the correlation functions we restrict to the Q range 0-2 GeV. The fits are, however, 
performed in the Q range 0-4 GeV. The dip around Q — 0.7 GeV and the following slope 
of the correlation function were treated as integral parts of the BE correlation function 
and as integral parts of the BEC effect. Therefore, all R2{Q) distributions are fitted with 
the parametrisation 

R2{Q) - 1 = Aii,e-««(1 + SdRQ) + 5^, (11) 

where Ar^ denotes the correlation strength, R is related to the source size, and the 
term with accommodates the dip around 0.6 < Q < 1.0 GeV. The auxiliary term 
6n accounts for small differences in the charged multiplicity of the signal and reference 
samples leading to a potential bias in the normalisation. 

The results of the correlation function, i?2(Q) — 1, are summarised in table 2 and shown 
in figure 4. In order to compare the measured correlation strength, Ar,^, between data 
and the BE32 model, the R and ea parameters were fixed to the values obtained from the 
model. The measured correlation strength in data was found to be slightly smaller than 
in the model. The significance of this difference is 1.8 standard deviations, but it does 
not include systematic errors. Fitting measurements of i?2(Q) must typically include an 
additional slope parameter and implies that it takes into account the uncertainty due 
to fragmentation. This leads to a large correlation with the parameter and means 
that it is no longer possible to extract meaningful results on e^, since it is not known 
whether to attribute the dip to the signal or the fragmentation of the reference sample. 
The R parameter will also be affected by this but to a lesser degree and finally the A^^ 
parameter is the most stable. No reliable way has been found to quantify the systematic 
errors on the R2{Q) measurements, and therefore only a qualitative agreement between 
BE32 and the semi-leptonic data can be demonstrated. 

When fitting 5i{Q) the normalization, S^, can no longer be described as a simple 
additive parameter, and is included in the fits as: 

5j{Q) = Aje-^'^{1 + eaRQ) + 5^{1 + -CMfO- (12) 
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Figure 7: The Ap(Q) (a,c) and D{Q) (b,d) distributions for like-sign and unlike-sign 
particle pairs, respectively. The BE32 Monte Carlo models including BEC between all 
particle pairs (BEA), pairs coming from the same W (BEI) are superimposed on the 
plots. 

The comparison of Si{Q) and R^iQ) results is hmited by the systematics of the R2{Q) 
measurements, as described above. Througout the paper, for reasons of clarity, the 
symbols A^j and A/ wil be used when referring to fits of the R2{Q) and Si{Q) distributions 
respectively. 

6.4 Statistical errors and correlations 

The values of the two-particle density distributions for different bins are statistically 
correlated, since in general a particle occurs in several pairs and because of non-Poissonian 
fluctuations in the overall particle multiplicity [21]. 

The covariance matrix elements for the histogrammed spectra are given by Vij = 
(hihj) — {hi){hj), where hi and hj are the contents of bins i and j. These covariances are 
propagated in the computation of the errors related to the distributions in equations (5)- 
(7). The use of pair weights does not pose any problems within this approach. The 
statistical errors shown in the figures in this paper were computed only from the diagonal 
elements of the covariance matrices. 

Because of the limited precision of V, the covariance matrix was treated carefully 
from a numerical point of view in the subsequent fits. The results and the covariance 
matrix were obtained from the same data and therefore correlated. This correlation is 
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Figure 8: The distribution for like-sign (a) and unlike-sign (b) particle pairs. The 

BE32 models including BEC between all particle pairs and pairs coming from the same W 
only are superimposed on the plots, together with the fit results using equation (12). The 
fit results are with fixed R and and are corrected for the bias mentioned in section 6.4. 

effectively enhanced by the substantial correlations in V and is a cause for biases in the 
fit results. By choosing suitable transformations of the fitting functions this bias can 
often be reduced. However, for the Si{Q) observable the value of A is already nearly 
unbiased when no transformation is applied. All these biases were found to be inversely 
proportional to the number of fitted events and the effect could therefore be estimated 
by comparing the fit results on large samples to the fit results on data-sized samples. 
The final biases on the parameters A and (5jv (see section 6.4) which were obtained for 
data-sized samples, were estimated to be 0.040 -1-0.031 ■ A and —0.0121, respectively. The 
biases on and R were found to be completely negligible. All the final fit results were 
corrected for these biases. 

The statistical errors were verified using resampling techniques and were found to be 
unbiased within a relative precision of 2%. 

7 Results 

Inter-W BE correlations were studied as function of the observables Ap{Q), D{Q) 
and 5i{Q) as defined in equations (5)-(7). These results are shown as function of Q 
in figures 7,8 and compared to predictions of the LUBOEI model. In the like-sign dis- 
tributions, an excess in data at low Q values can be observed. Numerical results were 
extracted from the 5i{Q) distribution, shown in figure 8. This choice was made since the 
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sample/pajrametcr 


A/ 


R{fm) 




Sn 




R free and fixed to BEA values 


BEA (±,±) 


1.50 ± 0.06 


0.72 ± 0.02 


-0.50 ± 0.03 


-0.010 ± 0.002 


116.4/96 


Data (±,±) 


1.42 ± 0.63 


1.14 ± 0.33 


-0.50 


-0.002 ± 0.020 


88.3/97 


BEA (+,-) 


0.30 ± 0.03 


0.41 ± 0.03 


-0.60 ± 0.09 


-0.010 ± 0.002 


110.8/96 


Data (+,-) 


0.43 ± 0.22 


0.45 ± 0.15 


-0.60 


0.000 ± 0.020 


93.2/97 


R and fixed to BEA values 


Data (±,±) 


0.82 ± 0.29 


0.72 


-0.50 


-0.005 ± 0.020 


89.9/98 


BEI (±,±) 


0.10 ± 0.05 


0.72 


-0.50 


-0.009 ± 0.004 


99.3/98 


BEO (±,±) 


0.02 ± 0.02 


0.72 


-0.50 


-0.010 ± 0.002 


98.0/98 


Data (+ ,-) 


0.40 ± 0.18 


0.41 


-0.60 


-0.001 ± 0.020 


93.2/98 


BEI (+ -) 


0.01 ± 0.03 


0.41 


-0.60 


-0.005 ± 0.003 


137.5/98 


BEO (+ ,-) 


-0.04 ± 0.02 


0.41 


-0.60 


-0.009 ± 0.002 


138.0/98 



Table 3: Fit results to like-sign and unlike-sign 5i{Q) with R free and it! fixed. Only 
statistical errors are shown. 

Si{Q) is a genuine inter-W correlation function having a clear physical interpretation. 
Other results are given for comparison. Note that none of the results are corrected for 
pion purity or secondary tracks, as was mentioned in section 3. 

The Sj{Q) distribution for like-sign pairs was fitted using equation (12). The fits to 
the BEA model were performed with all four fit parameters left free. The correlation 
strengths obtained from the BEI and BEO models were used to estimate a possible bias 
for the measurement which is treated as systematic error (see section 8). The result 
where both A/ and R were left free in the fit provides the most model-independent result. 
However, fixing the values of R and ea to the ones predicted by LUBOEI BE32 tuned 
to DELPHI inclusive Z data, makes model comparisons easier and contains the same 
statistical significance. 

The results on the data and models are summarised in table 3. The main result of 
this paper is then the observed correlation strength for like-sign pairs: 

Ai{±, ±) = 0.82 ± 0.29(stat) ± 0.17(syst), (13) 

with it! fixed to 0.72 fm. The first error is statistical and the second error is systematic. 
The evaluation of the systematic uncertainty of this result is discussed in section 8. The 
expectation from the BEA model yields 

Aj^^^(±, ±) = 1.50 ± 0.06(stat). (14) 

The correlations in unlike-sign pairs were also measured using the same procedure as 
for the like-sign pairs. The result is shown in figure 8 (b) and summarised in table 3. The 
distribution shows some enhancements at low Q when compared to the prediction from 
BEL Fitting with the same expression as for the like-sign pairs, but with R and e fixed 
to the unlike-sign BEA prediction, yields: 

Aj(+, -) = 0.40 ± 0.18(stat) ± 0.22(syst) (15) 

in agreement with the expectation from BEA: 



-) = 0.30 ± 0.03(stat). 



(16) 
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Figure 9: The 2-d Ax^ curves are shown for A/ and R. (a) shows the resuhs for hke-sign 
pairs, while in (b) the unhke-sign pair result is displayed. The three contours correspond 
to a Ax^ of 1, 4, and 9 respectively. The crosses show the prediction from LUBOEI 
BE32. The errors corresponding to the LUBOEI predictions are multiplied by a factor 5 
for clarity. 

The LUBOEI model prediction for unlike-sign pairs has not been tuned on Z events as 
it was the case for the like-sign pairs. Therefore this prediction should be treated more 
carefully and an interpretation will be made in the discussion section. 

The numerical results of the simultaneous fit of R and A/ are shown in table 3. Since 
the two parameters are strongly correlated the results are presented as Ax^ curves in 2-d 
plots. Figure 9 shows the fit results for both the like- and unlike-sign pairs. In these fits, 
the value of is still fixed since the data do not contain enough information about this 
parameter. The position of the dip in the correlation function is hence forced to scale 
like R. The systematic error is of the same size as the measurement with R fixed and 
not included in the fits shown in table 3 and figure 9. 

8 Systematic uncertainties 

The mo del- independent analysis of inter-W correlations performed in this paper avoids 
potential biases by the use of the event mixing and a direct data-to-data comparison. 
Therefore only few sources of systematic uncertainty remain to be considered; these are 
the subtraction of the Z — >■ qq background, the quality of the mixing of W final states 
from semi-leptonic events and the possible influence of colour reconnection. 

The background was subtracted using BE32 simulation where the correlation strength 
can be varied via the parameter PARJ(92). Although the standard BE32 tuning, using 
PARJ(92)=1.35, gives a good description of inclusive Z-events, it was found that an input 
parameter strength of PARJ(92)=0.9 gives a better description of Z-events having a clear 
four-jet topology. This can be seen in figure 10 where both tunings are compared with 
inclusive Z-decays and four-jet like events. The four-jet events were selected with the 
LUCLUS [17] clustering algorithm with ^4 > 4.0 GeV. 
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Figure 10: Top plots: The ratio of the two-particle densities of like-sign particle pairs 
for Z events in data and in Monte Carlo events using the BE32 model with parameters 
PARJ(92)=1.35 and PARJ(93)=0.34 GeV-\ a) for the inclusive sample, b) for a four-jet 
sample with ^34 > 4.0 GeV. The two bottom plots show the same comparison but with 
a different BE32 input parameter PARJ(92)=0.9. 

Both background samples were subtracted from the data and half of the absolute 
difference in the final result of A/, 0.075, was taken as a systematic uncertainty due to 
the lack of knowledge about BEC in qq(7) events. 

In addition, data with a lower purity than at the working point were used to estimate 
the correlations in the background itself. Four different purity bins were chosen and fits 
performed in each bin to obtain the dependency of A/ on the event purity. The mea- 
sured correlation strength depends non-linearly on the correlation strength parameter 
PARJ(92) used in the LUBOEI simulation of the background. A hnear interpolation be- 
tween the model predictions for the two parameter choices of the following form has there- 
fore been used: Model(6) = b- Model(PARJ(92)=1.35) + (!-&)■ Model(PARJ(92)=0.90). 
The extrapolation was found also to apply satisfactory outside the range < 6 < 1. 

In each purity bin and for several different background subtractions (as determined 
by different values of b) the data were then fitted in order to obtain the corresponding 
values of A/ and dAj/db. Only the uncorrelated errors were used in these fits, so that 
the results are independent except for the semi-leptonic events which are identical for 
each bin. The results of these fits are shown in table 4. The value of b was then varied 
to obtain the background subtraction, for which the A/ values are independent of the 
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Purity bin 


a 


b 


c 


d 


Purity 
A, {b = 0) 

dAi 
db 

Aj ih = -0.75) 


0.60 
0.39 ± 0.36 

-0.69 
0.91 ± 0.40 


0.77 
1.28 ± 0.35 

-0.48 
1.64 ± 0.36 


0.89 
0.74 ± 0.23 

-0.17 
0.87 ± 0.24 


0.97 
1.33 ± 0.33 

-0.08 
1.39 ± 0.33 



Table 4: Results for fits of individual bins in purity. Only the uncorrelated errors are 
shown. 

purity. The result of this variation gave the result: b = —0.75 ± 0.65, with = 4.1 for 3 
degrees of freedom. 

This value of b is consistent with the value used, b — 0.0 ± 0.5, and is an additional 
strong indication that the correlations in background 4-jet events are smaller than for 
the inclusive sample (6 = +1). The data are shown in figure 11, where the background 
with b = —0.75 is subtracted. The amount of subtracted background was varied by 10% 
(relative), which changed A/ by 0.007, which is neghgible. 

In the semi-leptonic events the background is a factor 2 smaller than in the fuUy- 
hadronic events. It was verified that the topologies of these backgrounds arc identical 
to those of the signal, i.e. mostly 2-fcrmion decays into hadrons. These events behave 
with respect to BEG in the same way as inclusive Z-events and therefore the BEG in 
these events are not expected to be significantly different than in the signal W-events. 
Finally, these background events do not suffer from the large extrapolation uncertainties 
of the 4-jet background in the fuUy-hadronic channel. The estimated systematic errors 
are shown in table 5. 

The selection of the data and the construction of the mixed reference sample may 
introduce distortions in the two-particle densities and therefore may lead to a non-zero 
value of A/, measured in Monte Garlo samples without inter-W BEG. The fragmentation 
models used at LEP do not give a perfect description of all the details of the soft frag- 
mentation and correlations. However, they constitute a reliable test for the magnitude 
of variation which can be expected for such effects. The largest absolute value of the 
measured Aj, 0.125, for these models (PYTHIA [17], HERWIG [22], ARIADNE [23]) 
was taken as a measure for the influence of selection procedures and mixing method on 
the measurement. The above method sums over many possible problems of the mixing. 
In order to verify that the sum is not small due to an accidental cancellation of large 
effects a weighting procedure was applied for several event variables and single particle 
distributions as described in section 4. BEI events were weighted so that the mixed and 
fully-hadronic events were in perfect agreement and the weighted events were refitted. 
The maximum shift in A/ was found to be 0.045 which is easily covered by the inclusive 
estimation. 

Detector effects are very small due to the analysis method. Any simple variation of 
detector performance was found to give zero shift in A/. Higher order effects are still 
possible, but these are also correlated with fragmentation properties and it was therefore 
assumed that they are covered by the previous estimate. 

A systematic uncertainty was attributed to the Golour Reconnection (GR) effect. This 
effect could have, in addition to inter-W BEG drastic consequences for the W mass 
measurement in the fully-hadronic channel [5,6]. As for BEG, it violates the assumption 
that the two produced W bosons decay independently of each other. 
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Figure 11: is plotted for the 4 different purity bins of table 4 (a-d). The background is 
subtracted using b — —0.75. The results are strongly correlated between the purity bins. 

Colour reconnection occurs when independent colour singlets interact strongly before 
hadron formation. In fuUy-hadronic W decays it recombines partons from different parton 
showers. After fragmentation, the resulting hadrons carry therefore a mixture of energy- 
momentum of both original W's. 

The CR effect has been modelled in various ways [5,6,24]. Only the extreme models [5], 
where reconnection occurs in all events, have been ruled out by the LEP experiments [9] ; 
however the absence of CR is also disfavoured by the same information. 

For this reason, three possible models of CR as implemented in JETSET, ARIADNE 
and HERWIG, were used to estimate their influence on this measurement. The maximum 
difference in A/ between the CR samples and their equivalent models without CR imple- 
mentation, 0.07, was found with the HERWIG implementation and was conservatively 
taken as systematic uncertainty due to the CR effect. 

Finally, half of the bias correction (0.033) described in section 6.4 was conservatively 
attributed as a systematic error due to fit biases. 

The total systematic uncertainty on the measured Aj value is the sum in quadrature 
of the contributions listed in table 5 for both like-sign and unlike-sign particle pairs. 
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syst source 


contribution to A7(±, ±) 


contribution to A7(+, — ) 


background BE model 


0.07 


0.02 


semi-leptonic Bg. 


0.01 


0.01 


cuts & mixing 


0.12 


0.04 


Colour Reconnection 


0.07 


0.22 


Bias Correction 


0.03 


0.03 


Total syst. 


0.17 


0.22 



Table 5: A breakdown of the systematic errors for the A/ measurement with fixed R, for 
like-sign and unlike-sign particle pairs. 

9 Discussion 

The result, A/ = 0.82 ± 0.29 ± 0.17, presented in section 7 shows an indication for 
BEC in like-sign pairs between two hadronically decaying W bosons at the level of 2.4 cr 
(standard deviations). The effect is 2.2 a below the prediction of BE32, assuming that the 
systematic uncertainties due to cuts and mixing are 100% correlated. The spatial scale, 
R, of the correlations is higher at the level of 1.3 cr, based on statistical uncertainties 
only. The QCD background to the WW signal is by itself of interest to study, and the 
data show that the BEC are weaker in these events than in inclusive hadronic Z events. 

The measurement where both R and Aj are left free and measured from the data 
contains the most model-independent information, while the result with fixed R is optimal 
for testing the excess in the direction predicted by BE32. 

The results for unlike-sign pairs are difficult to interpret due to a large model depen- 
dency and systematic errors. The unlike-sign pairs show a smaller excess at low Q values. 
The excess is at the 1.4 a level and in some agreement with the BE32 prediction, but a 
bit too large if the unlike-sign effect is scaled to follow the like-sign. The unhke-sign pairs 
are statistically correlated with the like-sign pairs at the level of 60%. Therefore it is not 
possible to rule out a large statistical component of these effects. 

The tuning of the BE32 model was verified using the semi-leptonic WW events. The 
overall agreement including the description of the dip around Q — 0.7 GeV was found to 
be excellent, except that the data show a bit less correlation than the prediction tuned 
on Z's. To this 1.8 o" statistical difference a systematic uncertainty which comes from 
extrapolating from Z data using the model has to be added. This uncertainty which is 
difficult to control is also the reason why the main result of this paper does not rely on 
BE32 simulation. 

Inside the LUBOEI models the input parameter PARJ(92) determines the strength of 
the correlations. However, due to non-linear effects in the implementation, there is no one- 
to-one correspondence between the generated strength and the measured strength, A^j /. 
Even with a fixed value of PARJ(92) the observed A^j does depend on the multiphcity of 
the event. The used value of PARJ(92)=1.35 leads to a measured value of A^^ = 0.77 ± 
0.02 for the semi-leptonic BE32 events while the BEA result is Aj = 1.50 ±0.06. The two 
BE32 results come from fits to R2 and Sj respectively. The systcmatics arc expected to be 
small on these numbers because they compare events which use the same fragmentation. 
The above observation can be combined with the information that the improvement in the 
statistical sensitivity to BE32 was smaller (9%) than the improvement (17%) expected 
from pure statistics. This indicates that the BEA not only adds correlations between 
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particles from different W's compared to BEI, but increases the strength of all correlations 
significantly. This enhancement scenario is not supported by the 4-jet Z-data where the 
BE32 model clearly overestimates the correlations when tuned on inclusive events. 

Both Ai^2,7 and R arc subject to additional corrections if their values are to be combined 
with other experiments. Af>,^j is diluted due to pair impurities to a level of about 70%, 
while R is only comparable to either other data or models. 

In conclusion this paper presents a model-independent measurement of inter- W BEC. 
The measurement does assume that intra- W BEC are the same in 2-jet and 4-jet WW 
decays. The direction in which the signal is looked for is motivated by the measured BEC 
in semi-leptonic WW events and implemented via the BE32 model tuned to Z data. The 
measurement is firmly based on the mixing of the hadronic final state of independent 
W bosons from semi-leptonic WW events. This technique yields small systematic uncer- 
tainties in spite of the small fraction of particle pairs from different W's. A weighting 
technique allowed the effective purity to be raised to about 20%. 

The treatment of statistical errors and bin correlations of the correlation function 
allows for a reliable specification of the statistical error of the quoted correlation strength. 
The remaining systematic uncertainties and their dependence on model parameters were 
carefully studied. 

10 Conclusion 

Overall, there is an indication of correlations in like-sign pairs from different W's 
with a significance of 2.4 standard deviations. The results are 2.2 standard deviations 
lower than the expectation of the BE32 model which was tuned to DELPHI Z data and 
verified on semi-leptonic WW events. The spatial scale of the correlations is larger but 
consistent with BE32. The results from the unlike-sign pairs are inconclusive, and their 
interpretation is model dependent. 

The prediction of the LUBOEI model is disfavoured not only by the behaviour of 
the WW data but also through the investigations of the background. It is clear that 
developments in the theoretical side of fragmentation models are needed before the BEC 
results presented in this paper can be fully understood. 

These final results of the DELPHI experiment have a hmited statistical accuracy. The 
precision would be significantly increased by a combination with the results of the LEP 
experiments which can be found in [20]. 
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